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Abstract 

The conditions for optimal reflection-free complex-absorbing potentials (CAPs) are discussed. It is 
shown that the CAPs as derived from the smooth-exterior-scaling transformation of the Hamiltonian, 
[J. Phys. B. 31, 1431 (1998)], serve as optimal reflection-free CAPs (RF-CAPs) in wave-packet 
propagation calculations of open systems. The initial wave packet, ^{t = 0) can be located in the 
interaction region (as in half collision experiments) where the CAPs have vanished or in the asymptote 
where VcAP 7^ 0. As we show the optimal CAPs can he introduced also in the region where the physical 
potential has not vanished. The un-avoided reflections due to the use of a flnite number of grid points 
(or basis functions) are discussed. A simple way to reduce the "edge-grid" reflection effect is described. 
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1. INTRODUCTION 



There is an extensive use in wave-packet (WP) propagation calculations in complex absorbing 
potentials (CAPs). The use of CAPs in propagation of WP calculations is usually for half 
collision experiments where the initial wave-packet is localized in the interaction region where 
the CAPs are vanished. The role of the CAPs is to avoid the reflection from the edge of the grid 
as obtained in the numerical propagation calculations. Often CAPs are referred to as optical 
potentials. The CAPs are used in very different fields of physics, chemistry and technology. See 
for example: calculations of resonances for CAPs in a nuclear physics problem^; deriving new 
expressions that simplify the numerical calculations of state-to-state transitions probabilities 
for reactive scattering collisions (for time independent Hamiltonians see Ref. and for time 
dependent ones see Ref. jsj) ; calculations of complex molecular potential energy surfaces by 



CAPsj^; and molecular electronic studies where the CAP serves to absorb charge reaching the 
electrodes^]. Beside the use of CAPs in the numerical calculations an effort has been taken in 

developing different type of CAPs. See for example Ref. B fl y, y, 0, fl Q where recently 

new type of CAPs were developed. For a most recent review on CAPs see Ref. jisj] . 

In 1998 we have derived CAPs by applying the smooth-exterior-scaling transformations 
(SES) to the HamiltonianQ. Here we study the conjecture that the use of exterior-scaling or 
SES similarity transformations, produce reflection free CAPs (RF-CAPs) for the WP propaga- 
tion calculations. As we will show here within the finite basis-set or finite grid approximations 
the CAPs are not reflection free ones. However, it is possible to show that for a given finite 
basis/grid method a quantity criteria for the strength of the numerical reflections can be de- 
rived. It is important to mention that about the same time Riss and MayerQ] obtained CAPs, 
which under specific conditions are similar to the SES-CAPs, by taking another approach (so 
called Transformative-CAP). Only when the CAP is introduced in a region where the potential 
energy has been vanished, the Transformative-CAP derived by Riss and Meyer is equal 
to the SES-CAP that has been derived analytically without any approximations by usQ] (in 
such a case the SES-CAP and the Transformative-CAP are identical although they look slightly 
different). Our main motivation for deriving SES-CAPs was to simplify the calculations of res- 
onances positions and widths. However, this SES-CAP has been used also to avoid the artificial 
reflections from the edge of the grid in wave-packet (WP) propagation calculations jl^. One 
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may wonder, what is the need for the Transformative- CAP s or the SES-CAPs since the reflec- 
tions can be taken as small as one wishes by introducing the CAP in the domain where the 
physical potential is zero and by making the CAP (any CAP) soft and long enough [itIi. The 
answer to that question is that it is most desired to avoid the use of long ranged CAPs which 
require large number of basis functions or large number of grid points in heavy duty numerical 
calculations. For example, in propagation calculations of many electron molecular systems it 
is hard to avoid the introduction of the CAP in the domain where the physical long range 
potential is not zero. 

Here we want to discuss two type of questions. The first type are mathematical-physical 
questions (i.e., theoretical questions in the sense that we assume that complete basis sets are 
used). Such as, 

- What are the properties of reflection free CAPs ? (as we show here, there are two conditions 
that should be satisfied). 

- Can we introduce the RF-CAP in the domain where the physical potential is not zero and 
the propagated wave packet does not consists of out going waves only? (as we show the answer 
is yes). 

- Can the initial state be exponentially localized in the interaction region, as required in 
half-collision experiments, where the CAP vanishes? (the answer is yes). 

- Can the initial state be localized in the domain where Vcap 7^ ? (the answer is yes, 
provided the smooth-exterior scaling transformation is applied to the initial state). 

The second type are practical questions. 

- Are indeed the RF-CAPs reflection free in the numerical calculations where finite number 
of grid-points or finite of basis functions are used ? (the answer is no since in spite of the 
complete absorbing of the fast moving components of the wavepacket still there is an edge-grid 
reflection effect which is associated with the slow moving components of the WP). 

- Can we minimize the reflections which result from the use of finite sized basis/grid methods 
and how ? (the answer is yes, by methods explained in the paper). 

- Can we apply the RF-CAPs to many electron problems ? (the answer is yes provided the 
electronic repulsion terms, l/\fi — fj\ are modified. This requirement can be avoided when the 
ionized electrons are not correlated). 
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2. WHAT ARE THE IDEAL REFLECTION-FREE CAPS ? 



First, we should describe the numerical problem we want to solve by introducing a CAP 
into the Hamiltonian. Using the hermitian quantum mechanics the propagated wave packet is 
given by, 



In the numerical calculations the propagated wave packet is $num(^) 7^ ^exact{t)- We are 
looking for numerical methods for which. 



where e determines the requirement accuracy from the numerical results. Since in the numerical 
calculations only finite number of grid points or finite number of basis functions are used, the 
available spatial space is not from r = to r = oo but up to r = L. Therefore, accurate 
results are obtained as long as ^exactit) vanishes at r > L. By increasing the number of the 
grid-points or by increasing the number of the basis functions we increase the value of L. The 
initial state, $0) is a square integrable function. In half collision experiments (such in photo- 
dissociative or auto-ionization reactions) the initial WP is localized in the interaction region 
where l^numit = 0) — ^exactit = 0) | < €. However, as time passes the wave-packet spreads and 
only during a given period of time r, the numerical calculations satisfy the accuracy condition 
stated above. It is important to realize that the value of r is determined by the time it takes 
for the tail of the wave-packet to reach the edge of spatial space (i.e., r=L). In order to obtain 
^numit) within the desired accuracy, one should increase the number of the used grid/basis 
points/functions and thereby increase the value of L. The role of the CAP is to enable one to 
obtain accurate numerical results in the limited available spatial space, r < rcAP < L, without 
the need to increase the number of grid/basis points /functions. Namely, 



^ exactif^^ 



(1) 



exact 



(2) 



(3) 



where due to the use of the finite grid/basis-set numerical methods. 



$CAp(r >L,t) = 0. 



(4) 



The CAP is defined such that. 



VcAP = as r < rcAP < L 



(5) 
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and, 

l^cApit) - ^exact{t)\ < t 171 the regioTi where Vqap = 0. 
A common requirement is that, 

$0 = in the region where Vcap 7^ 



(6) 



(7) 



As a matter of fact the last condition is too strong and it is possible to satisfy EqlHl also when 
the initial state is localized in the region where the CAP gets non zero values. This extension 
will be discussed later. 

Short range CAP (the Saxson-Wood potential) has been used about two decades ago in 



molecular wave- pac ket calcu^ 
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18 L A CAP which has been used often in the literature 
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is Vo = for X < xo and 



Vo = — '^A(x — Xo)" where n = 1,2, ...8 for x > 0. For large values of n these CAPs are very 
similar to the purely imaginary step type potential that has been shown above to provide a 
strong reflection. Regarding the reflections due to the introducing of abrupt complex potentials 
one might be aware to the fact that there are examples (see the review in Ref . and references 
therein) of discontinuous potentials that are constructed to avoid reflection, and absorb totally, 
at single incident energies, or in certain momentum intervals, or at a discrete set of energies. Of 
course they cause reflections at other energies. The CAPs that we are looking for are different 
ones. They are energy independent RF-CAPs, and in principle can be chosen to be universal 
ones (i.e., problem independent). 

As we will show here it is unlikely to have a universal (i.e., problem independent) CAP for 
which both EqlHand EqlHlare satisfied. Therefore, let us first discuss the possibility to satisfy 
Eqiniwhen the condition given by Eq|3]is replaced by a weaker numerical condition: ^cAp{t) 
is a square integrable function at any given time, which decays to zero much faster than the 
exact solution. Such that within a given time interval, 

^cApir = L{T),t<T)<e, (8) 



where the value of e is determined from the desired accuracy of the numerical calculations. 
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3. THE SES-TRANSFORMATIONS AND THE CONDITIONS FOR OPTIMAL 
REFLECTION-FREE CAPS. 



The idea of introducing RF-CAPs by using the exterior scahng or SES methods is clear: 
the Hamiltonian remains as it is inside the inner region, where the coordinates stay on the real 
axis. However, it has been shown by Simon that upon the exterior scaling transformation, 

r Text (9) 

where inside the inner unsealed region, 

Text = r when r < tcap (10) 

and in the external-scaled region, 

rext = rcAP + {r - rcAp)e''^ when r > tqap (11) 

the eigenf unctions are not equal to eigenf unctions of the unsealed (i.e., hermitian) problem 
inside the unsealed region j35|. For example, for a free particle Hamiltonian the contin- 
uum eigenf unctions inside the inner unsealed region are given by. Am exp{—ik exp{—i6)r) + 
Ao^i exp(+iA;exp(— 'j^^)r). Since the propagated WP can be described as a linear combination 
of the eigenf unctions of the complex scaled (or exterior scaled) Hamiltonian, it is not clear 
at all that in this case EqlH] is satisfied (here we consider the exterior scaled Hamiltonian as 
H + VcAp)- This result is very confusing since from numerical propagation calculations we 
know that inside the inner unsealed region in space, the propagated WP is exactly as obtained 
without the use of exterior-scaling. As we will show below the validity EqE] can be easily 
explained by association the SES approach with the use of similarity transformation operators 
as developed in Ref . [l^ Isfij] . Using the SES approach, 

r ^ rsES = Fe{r), (12) 

where the path in the complex coordinate space is chosen such that, 

\Fg{r) — r\ < e when r < tqap (13) 



and. 



Mi ^ as r^oo. (14) 
r 
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The SES transformations clearly show that EqEl can be satisfied to any desired accuracy. If 
the SES transformation is represented by the similarity operator, S, than the propagated WP 
within the framework of the SES approach is given by S'^exactit) which is equal to "^exactit) 
inside the inner region (see Eq]13p where S* ~ 1. 

Let us discuss now the validity of Eq|H] Following Simon's proof for the exterior 
scaled potential and following Moiseyev and Hirschfelder's proof for general complex scaled 
transformationsjs^l (including the SES transformations), the complex scaled resonances func- 
tions are square integrable but the continuum eigenfunctions are not square integrable functions. 
They are associated with complex eigenvalues, E(..j.t{continuum) = k1^^/2 = {k expl—iO)"^ /2. 
Such that fcezi ■ "^ext in the exterior region is equal to the same value as obtained in hermitian 
quantum mechanics, i.e., k ■ r (note that a very different result is obtained in the inner region as 
discussed in the previous paragraph). Therefore, the asymptote of the continuum wavefunctions 
as obtained after the application of the exterior or SES transformations, remain as obtained 
within the framework of the conventional (i.e., hermitian) QM approach. Upon complex scaling 
kext is rotated into the lower-half complex k-plane to avoid the exponentially divergence of the 
complex scaled incoming waves associated with real and positive values for the wave vector, i.e., 
exp{—ikrext) = exp{—ikcos{d)r) exp{+ksm{6)r) —>■ oo, as r oo. Therefore, it is no obvious 
weather a square integrable WP such as, ^exact 

(t) = J^^C{k,t)exp{tkr)dk as r — > cxD, 
remains square integrable when r — >■ Vext or r — > tses- It has been proven by Moiseyev and 
Katrieljs^ that for sufficiently small values of 9, i.e., 9 < 9c, the eigenfunctions of a complex 
scaled Hamiltonian which are associated with the bound states are square integrable. The 
value of 9c depends on the shape of the potentialfs^. Let us assume that the wave-packet is a 
Gaussian, exp(— ar^). It is clear that exp{—ar'^^f- g^g) remains square integrable provided that, 
9 < 9c = tt/2. When the wave-packet is more localized, for example is described as exp(— ar^), 
then 9c = tt/N. Since Gaussians form an over-complete basis set, one might expect that any 
square integrable function (which can be expanded in term of the Gaussian basis set) remains 
square integrable after applying the complex scaling or the SES transformation. 

It is easy to prove that the wave-packet ^exact{t) decays exponentially to zero at any given 
time, provided it is a square integrable function at t=0. A proof which holds also for complex 
scaled non- hermitian Hamiltonians is as follow: ^exact(t + dt) = exp{—iHdt /h)^exact(t) . For 
sufficiently small value of dt, exp{—iHdt/h) = X]n=o(^0~"'^(~"'^'^'^/^)"(-^)" is ^ converged series 
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(provided ImH < 0). If ^exact{t' < t) is a square integrable function than {HY'^exactit) 
is square integrable as well (the second derivative of a square integrable function is square 
integrable and the product of a square integrable function and a confined (complex scaled) 
potential is also a square integrable function). 

Let us summarize the facts we know by now: (1) when the initial wave-packet (WP) is 
square integrable the time propagated WP is square integrable as well; (2) the complex scaled 
square integrable WP remains square integrable; (3) the complex scaled incoming waves diverge 
exponentially whereas the outgoing waves exponentially decay to zero; (4) in the absence of 
a source of particles in infinite large distance from the studied system, the asymptote of the 
propagated WP consists of out going waves only (as in half collision experiments). From (l)-(3) 
it is clear that for the most general case the square integrable WP, 

^exact{r >L,t)= / dk{D{k, 1)6"'^' + C{k, t)e+^'=") ^ as r ^ oo, (15) 
Jo 

remains square integrable, 

/■oo 

^CAp{r>L,t) = / dk{D{k,t)e-'^^"''^^^''e^^'^^^^'' + C{k,t)e^'^^''''^^^''e (16) 
-H^ as r oo, (17) 

although each one of the components of the complex scaled incoming waves exponentially 
diverge. When the condition (4) is not satisfied this fact (i.e., interference of exponentially 
diverged incoming waves results in a square integrable function) may introduce some numerical 
difficulties in the propagation calculations. For overcoming these type of numerical difficulties 
when long ranged potentials are used see the second reference in Ref. j^]. 
When condition (4) is satisfied (as in all half collisions experiments) then, 

/■oo 

^exact{r >L,t)= / dkC{k, t)e+*^ _ as r oo, (18) 
Jo 

and it is easy to see that $cap decays faster since, 

POO 

^CAp{r>L,t)= c^A;C(fc,^)e+''='=°^(^)"e-'=^'"(^)^ (19) 







The fact that within the interval of tcap < r < L, the propagated WP, ^cApij), decays faster 
than ^excat{t), is the main motivation behind the use of the exterior scaling, smooth-exterior 
complex scaling methods in the numerical propagation calculations. 
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4. A QUANTITY CRITERIA FOR THE MEASUREMENT OF THE STRENGTH 
OF THE NUMERICAL REFLECTIONS 



From Eq]19la quantity criteria for the strength of the numerical reflections from the edge of 
the grid is obtained, 

|C(A;,t)e-'=""(^)^| < e. (20) 
As an upper limit for the accuracy of the calculations one gets that, 

\C{k,t)e^''^\ < e. (21) 

At t=0 the initial wave-packet gets exponentially small values at r > L and therefore we 
can consider it where, C{k, 0) = 0. As time passes the wavelet with the largest value 

of k (associated with a large velocity) is the first to reach the edge of the grid. As one can see 
from Eql^the fast moving components of the wavepacket are entirely absorbed at r = L, due 
the use of the complex absorbing boundary conditions which were introduced by the use of the 
exterior scaling or the SES transformations. For the components of the wavepacket associated 
with small values of k, the requirement of exp{—kL) ~ is satisfied by increasing the value of 
L. The propagation calculations using SES transformations, within the framework of the finite 
basis-set /grid approximations, are accurate as long as \C{k,t)\ gets sufficiently small values. 
This explains why L in EqISl is a function of time and why the duration of the propagation 
calculations can not exceed a given period of time, T, when L is held fixed in the propagation 
calculations. 

As an illustrative example we carried out wavepacket propagation calculations for a one- 
dimensional Gaussian, \l/(x, t = 0) = (l/Svr)^/^ exp(— + ipo^;)) which is localized at a 
potential well embedded in between two identical potential barriers. This potential, V{x) = 
(0.5x^ — 0.8) exp(— O.la;^), has been used before as a test problem for new methods developed in 
non-hermitian quantum mechanics (see for example and references therein). In Fig. 1 the 
results obtained from two type of propagation calculations are presented. The long dashed line 
stands for the numerically exact calculations of "^exactixjt = 60), using 5-order split operator 
with —1000 < X < +1000. The full sohd line is \I/(-Fe=o.5rad(a;), t = 60) where Fe(x) is a smooth 
exterior scaling function. Such that, Fe ~ a; when |x| < xcap = 90, whereas Fq = xexp{i9) 
when |x| > xcap- For |x| < xcap '^exact(yX,t = 60) ~ '^{Fg{x),t = 60). However, it is clearly 
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FIG. 1: The numerical exact propagated wavepacket (long dashed line), ^exact{x,t = 60), and the 
corresponding wavepacket (denoted by NH-QM) which is defined as '^{Fg=Q,5rad{x),t = 60). The 
smooth-exterior-scaling contour is defined as, Fg ^ x when \x\ < xcap = 90, whereas Fg = xexp{i9) 
when |x| > xcAP- The initial wavepacket is given by, ^(x,t = 0) = (l/57r)^/''exp(-xVlO + ipoa;), 
where po = 1. 

shown that unlike the exact wavepacket which oscillates, the smooth exterior scaled wavepacket 
(labelled in Fig.l by non-hermitian quantum mechanics (NH-QM) approach) decays to zero as 
X is rotated into the complex coordinate plane around x = xqap = 90. 

Following our analysis the propagated wavepacket decays to zero when the contour x is 
smooth exterior scaled (rotated) into the complex coordinate space only within the time interval 
t < T. The results presented (denoted by NH-QM) in Fig. 2 were obtained from numerical 
calculations where —100 < x < +100 (i.e., the box size is L = 200). It is clearly shown that 
until t < 30 the complex scaled wavepacket is practically equal to zero at the edge of the grid 
(i.e., at x=L/2). As time exceeds the value of t = T = 30 the complex scaled wavepacket is 
reflected from the edge of the grid. 

5. HOW TO REDUCE THE NUMERICAL REFLECTIONS OF THE SLOW MOV- 
ING COMPONENTS OF THE WAVEPACKET FROM THE EDGE OF THE GRID 

Let us propose two different possibilities: 

(a) Accelerate the slow moving components of the wavepacket by inducing an external dc-field. 
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L=200 




NH-QM/dc 



FIG. 2: The non-hermitian propagated wavepacket which is constructed from 400 Fourier basis 
functions (long dashed hne, denoted by NH-QM) as function of time at x = L/2 (edge of the grid). 
The propagated wavepacket is defined as ^'(-Fe=o.5rad(3^ = 100), f). Such that, Fg ^ x when \x\ < 
xcAP = 90, whereas Fq = xexp{i9) when |x| > xcap- The initial wavepacket is given by, ^'(x,i = 
0) = (l/57r)^/^ exp(— x^/10 + ipox), where pq = 0. The reflections from the edge of the grid as time 
passes are obtained when t > T = 30. The full line stands for the results obtained when a dc-field has 
been introduced close to the edge of the grid, x^f. = 95 and £dc = 2. 

Eql^ indicates clearly that the numerical reflections from the edge of the grid, are associ- 
ated with slow moving wavelengths. As discussed above the fast moving components of the 
wavepacket are entirely absorbed at r = L, due the use of the SES-CAPs. A possible solution 
to this problem is by adding a static field close to the edge of the grid, in order to accelerate the 
slow moving wavelengths which are completely absorbed by the SES-CAPs. The static field is 
turned on only at the edge of the grid where r > and is given by, 

Vext-dc{r > Tdc) = ~ ^dc)- (22) 

The value of r^c and Sdc can be optimized to minimize the effect of the dc-field on the WP 
propagation. An estimate of the error introduced by adding the Vext-dc potential term to the 
Hamiltonian can be obtained from the imaginary parts of the bound states calculated for the 
Hamiltonian which is taken as, H + Vext-dc + Vcap- The evaluation of Vcap, from the SES 
transformations will be described in Section 6. 

Let us return to our illustrative numerical example. The smooth exterior scaled wavepacket 
has been calculated as described above when |a;| < 100 (i.e., the box size is L = 200 ) when 
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FIG. 3: The numerical propagated wavepackets at t=250 as obtained when 400 Fourier functions 
were used as a basis set. The reflections from the edge of the grid (which appear at t > 30 as shown 
in Fig. 2) are avoided as a dc field is added close to the edge of the grid. 

Vext-dc dc-potential term was added into the Hamiltonian. The results presented in Fig. 2 (the 
full solid line denoted by NH-QM/dc) clearly show the strong suppression of the reflections 
from the edge of the grid as the dc field has been introduced into the Hamiltonian. In Fig. 3 
we present the numerical results obtained at t = 250. It is clearly shown that close to the edge 
of the grid the dc field inhibits the artificial reflections as appeared in the NH-QM calculations. 
Note that as time passes the reflection leads to the distortion of the wavepacket also at regions 
which are quite far from the edge of the grid, as it is shown in Fig. 4. However, as one can 
see from the results presented in Fig. 4 the introducing of the static fleld reduces this artiflcial 
edge-grid reflection effect. 

(b) Imposing of out- going boundary conditions 

The numerical edge-grid reflection effect can be reduced by imposing out-going boundary 
conditions (complex scaled ones in our case). It is simple to implement that approach when 
grid methods are used, 

$CAp{t + At)=lJ$cAp{t), (23) 

where U(t + At ^ t) is the time evolution NxN matrix and the jth component of the vector 
^CApit) is the value of the propagated WP at fj;j = 1, 2, ..A^ grid point. The grids points are 
ordered such that |ri| < |r2| < ...|r7v_i| < |rAr|. The notation of CAP stands for the use of 
the exterior or the SES transformations. The out-going boundary conditions are imposed by 



12 




FIG. 4: The propagated wavepackets at t=250 as obtained when 400 Fourier functions were used 
as a basis set, in a comparison with the numericahy exact solution. The propagated wavepacket as 
obtained when both the smooth exterior scaUng transformation and the dc field were implemented 
into the numerical calculations, is in a very good agreement with the numerical exact solution. 

replacing the N-th row of the time evolution matrix UN,i', ^ = 1, 2, by exp[i exp{i6)k{fN — 
'i^N~i)]UN-i,i- Here we use the fact that the r^r, -r/v-i, 'nv-2 grid points are in the scaled region 
where r — fexp{i6). The wave vector k is determined from the previous time step calculations 
and from the tn-i and the r/v_2 grid points. That is, 

exp[^exp(^^)fc(r-W-l - fN-2)] = ^"^^"^i?"'' ■ (24) 

It should be stressed that it is not always true that at a given time the tail of the WP is 
constructed of a single out-going wave component. However, this kind of an approximation has 
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been found useful in WP propagation calculations of various physical problems 

One should assure that the real part of the wave vectors get positive values only. In the one- 
dimensional case where equally spaced grid points are used, the application of that approach is 
straightforward. In such a case the modified last row of the time evolution matrix is given by, 

f/ (t + At^t) = UM-iAt + ^ ^)ir^^P^^v (25) 

^CAp[rN-2,t) 

Similarly, one can modify all j > jc rows and not only the last one. The assumption is that 
the vectors, rj^j^ are all embedded in the asymptote region of the propagated wavepacket. 
This method (applicable to 3D problems as well) to reduce the edge-grid reflection effect is an 
extension/variation of Hadley's original work, where the transparent boundary condition for 
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FIG. 5: The propagated wavepacket as obtained when the time evolution operator has been modified 
as explained in the text (see Ea l25() . The model Hamiltonian and the initial state are as described 
in the caption of Fig. 1. The long dashed line stands for the results obtained from numerically exact 
propagation calculations. 

beam propagation method was developed 2^ . This method does not require the use of CAPs. 
However, we believe that the use of the exterior or the SES transformation together with the 
transparent boundary condition, should minimize the numerical reflections from the edge of 
the grid. 

The possibly to impose outgoing boundary condition by modifying the time evolution oper- 
ator as shown in Eq O^ is illustrated here by applying it to our test-case problem. The results 
presented in Fig. 5 clearly show that by using the method introduced above a similar absorbing 
boundary condition effect - as achieved when the RF-CAPS are added to the Hamiltonian - 
is obtained. A comparison between the results presented in Fig. 5 and Fig. 1 shows that the 
RF-CAP (Fig. 1) provides better results than those obtained by imposing out going boundary 
condition on the propagated wavepacket (Fig. 5). However, it might be expected that similarly 
to the effect of the dc-field described above, the combination of the two approaches would avoid 
the reflections which are obtained after long time propagation. This study is out of the scope 
of the present study and requires a further investigation. 
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6. THE OPTIMAL RF-CAPS FROM THE SES TRANSFORMATIONS. 



In order to complete the representation of optimal RF-CAPs we should show that under well 
defined specific conditions the use of the SES transformation, is equivalent to the inclusion of 
a CAP which gets non-zero values only in the edge of the grid. This SES-CAP is an non-local 
operator, since it includes terms with the momentum and kinetic energy operators. For the 
sake of a coherent representation of the subject we briefly describe how the SES-CAPs are 
obtained. 

The SES transformed time dependent Schrodinger equation can be rewritten as, 

d 

UcAP^CApit) = i — ^cAp{t) (26) 

where, 

HcAP = SHS-' (27) 



$CAP = S^e.act{t) = $e.act (^^9 (r) , t) (28) 

We have proved before that the SES transformation is equivalent to the including of a 
non-local energy independent, universal (i.e., problem independent) CApjl3|. 

HcAP = H + AV + Vrf-cap, (29) 

where the correction term to the physical potential is given by, 

AV = V{Fg{r))-V{r), (30) 



and the non-local energy independent, universal CAP has been proved to be equal to 



Vrf-cap = Vo{r, 9) + V{r, 9)— + V^ir, 9)— (31) 

or or'^ 

The complex functions V^- ; j = 0, 1, 2 are vanished in the inner region where r < tcap- They are 
inverse proportional to the mass (reduced mass) of the particle (s) which is (are) absorbed by 
'he Vrf-cap and are defined as ( note that below we use the notation Fg^\r) = dFg{r)/dr^) 
3: 



Voir, 9) 
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.(3)..^ 5(^f W)^ 



2Fi'^(r) 



(32) 



V,{r,e) = ^\^ (33) 
M(F,(')(r))3 

y2{r,e) = ^(l-{F!,'\r)r') (34) 



2M 

The initial state is defined as, 

(^oitrans formed) = $(Fe(r),t = 0) (35) 

When the initial state is localized in the interaction region where Vrf-cap = then 
$(Fe(r),t = 0) = $(r, t = 0). In the case that the physical potential is a short range po- 
tential it is quite obvious that the contour of integration Fg{r) can be chosen to yield AV = 
everywhere at any point in the entire space. In the case of long range potential the situ- 
ation is more complicated 13]. In such a case the AV term in Eq|^ can not be neglected 
and the SES RF-CAP is equal to AV + Vrf-cap and seems to be problem dependent. How- 
ever, for neutral molecules if rcAP gets a sufficient large value such that the ionized elec- 
trons are in hydrogenic like orbitals then AV + Vcap can be replaced by a universal poten- 
tial term 1/r — 1/Fg{r) + Vrf-cap 3|- This approach holds also for many electron systems 
where we assume that the ionized electrons are not correlated as they get far away from the 
atom/molecule/QD. In such a case we do not need to replace the two electron repulsion terms 
|fj - fj\-^ by \Fg{ri) - Fg{rj)\-^. 

Before concluding let us return to our illustrative numerical example. Using 400 Fourier 
basis functions (with the box-size, L=200) we obtained matrix representations of the hermitian 
Hamiltonian, Hh-qm = — O.Srf^/dx^ + (0.5a;^ — 0.8) exp(— O.lx^), and also of the non-hermitian 
one, Hn h- om = Hh-qm + Vrf-cap- The parameters for the function Fg{x) as defined in 
Ref. [ij, |l6| are 9 = 0.5rad, A = 0.9 and xcap = 90. The two 400x400 matrices where 
diagonalized. The eigenvalues and eigenvectors of the Hermitian matrix are correspondingly 
given by, E^^^'^^real) and Cf'^'^^. Similarly, E^^~^^\complex) and (j^^-Q'^'^ are associated 
with the eigenvalues and eigenvectors of the non-hermitian (complex and symmetric) matrix. 
The propagated wavepacket within the framework of H-QM is given by, 

400 ^ 

'^H-QM{x,t) = 2_^e ^ 2^C^,j exp(z— ^) (36) 

j=l n 
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t=60 




X 



FIG. 6: The propagated wavepackets at t=60 as obtained from conventional and non-hermitian QM 
calculations, in a comparison with the numerically exact solution (denoted by a long dashed line). The 
propagated wavepackets denoted by H-QM and NH-QM correspondingly were constructed from the 
eigenfunctions and eigenvalues of the hermitian and non-hermitian Hamiltonians when 400 Fourier 
basis functions were used as a basis set. The NH-QM results obtained when the smooth-exterior- 
scaling transformation was introduced ( by adding Vrf-cap into the Hermitian hamiltonian) , are in 
a complete agreement with the exact solution when |x| < xqap = 90. 

The propagated wavepacket within the framework of the NH-QM approach is given by, 

400 „ 

'^NH-QM{x,t) = 2_^e J Z^'^nj exp(z— ^) (37) 

j=l n 

The resuhs presented in Fig. 6 clearly show that while the reflections from the edges of the 
grid appeared in the propagation calculations within the framework of the conventional QM, 
they do not show up in the NH-QM calculations. The reflections appeared in the conventional 
quantum mechanical calculations due to the use of the eigenfunctions, which were obtained 
within the framework of the box-quantization approximation, as a basis set. A quasi-discrete 
continuum rather than a continuous continuum has been used in the propagation calcula- 
tions. As one can see from the results presented in Fig. 6 the use of the RF-CAP provides, 
^7VH-QA/(-90 < X < +90, t = 60) = '^exact{-^^ < X < +90, t = 60) (within more than 6 
digits of accuracy). In spite of the fact that in the two calculations we have used the same 
basis functions and the same number of them the NH-QM calculations provided an accurate 
propagated wavepacket while the conventional calculations are far from convergence. 
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7. CONCLUDING REMARKS 

rn 

We can summarize it by saying that for the CAPs derived from the SES transformations (14] : 
(1) the propagated WP decays faster to zero than the exact solution and therefore at any given 
time we can use a smaller grid/basis in the numerical calculations when the SES-CAPs are 
introduced into the numerical calculations; (2) the SES-CAPs can be introduced also in the 
region where the interaction potential is active (provided the edge of the grid is in the region 
where the exact WP has outgoing wave components only); (3) the use of SES-CAPs enables 
one to introduce the CAPs also in the region where the initial WP does not get zero values; 
(4) the duration of the WP calculations which provide accurate results (avoiding the numerical 
reflections from the edge of the grid) can be easily estimated (see Fig. 3); (5) it is possible to 
reduce the reflections of the slow moving components of the wavepacket, either by introducing a 
dc-field in the edge of the grid or by imposing out-going boundary conditions on the propagated 
WP. (6) The SES-CAPs are indeed the optimal reflection free caps, RF-CAPS, for wavepacket 
propagation calculations. 
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